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1 Introduction 

We consider the problem of waves of the free surface of an inviscid and incompressible fluid 
which satisfies the Euler equation. It is well known that this problem can be formulated as 
an Hamiltonian system (Zakharov, 1968) with canonical conjugate variables (77 (x), -0(2;)). 
Here r){x) is the free surface displacement above the point x of the (planar horizontal) 
bottom, and ip{x) is the velocity potential taken at the surface ip{x) — ^{x, r){x)) . There 
are infinitely many degrees of freedom, labeled by the parameter x . 
Several important numerical and analytical studies have been performed in connection 
with the question of existence and stability of travelling (or permanent form) waves. 
They are solutions of the form rjlx, t) — g{x — c.t) for some function g and constant c. 
We refer in particular to the works of (Levi-Civita, 1925), (Crapper, 1957), (Longuet- 
Higgins, 1978), (McLean, 1980), (Chen & Saffman, 1980), (Mackay & Saffman, 1986), 
(Craig & Sulem, 1993), (Craig & Worfolk, 1995), (Craig, 1996), (Bridges & Bias, 1996) 
and (Debiane & Kharif, 1996). 

Prom now on, we restrict ourself to the case where the variable x e R, and the solutions 
T],'^ are required to satisfy periodic boundary conditions, in x (in [0, 27r] for instance). 
Of particular interest for searching the solutions of an Hamiltonian system, the Birkhoff 
reduction to normal forms (Birkhoff, 1927) is known to be a powerful method of analysis. 
It even gives rise, for some cases of the water waves problem, to integrable reduced 
Hamiltonians as it was shown in (Dyachenko & Zakharov, 1994). However in (Craig & 
Worfolk, 1995), the authors shown that contrary to the conjecture stated in the preceding 
paper, the fifth order Birkhoff form of the full Zakharov Hamiltonian is no more integrable 
due to the presence of resonant terms. It is well known that such terms will prevent the 
Birkhoff iterative process to converge (problem of small divisors): see for instance (Siegel 
& Moser, 1971). 

On the other hand, the KAM theory, (Kolmogorov, 1954), (Arnold, 1963) and (Moser, 
1963) indicates that under some conditions, an Hamiltonian system treated as a quasi- 
integrablc one exhibits a great number of quasi-periodic solutions on invariant tori, 
provided that the integrable part of the Hamiltonian system is non degenerate, and that 
the perturbation is small enough. By quasi-integrable we mean that the Hamiltonian can 
be written as an integrable one, added with a small quantities in a sense to be precised. 
We also need that some parameters (or initial conditions) verify a non-resonance condition 
of Diophantine type: i.e. they must be in some Cantor set. 

In this paper we use the KAM method to prove the existence of a solution of a simplified, 
non-integrable model of surface gravity waves in deep water, introduced by (Zufiria, 
1987). As usual in KAM theory, we build a canonical transformation that takes the 
initial Hamiltonian into an integrable one, at least locally in the phase space. In this 
way it is possible to write a trajectory of the initial Hamiltonian system for each initial 
non-resonant condition. It is easy to show, by inspection, that such trajectories are not 
travelling waves. It will be clear that the method can be used for any other truncated 
(finite dimensional) form of the Zakharov Hamiltonian. 

This paper is organized as follows: 

In section 2, we describe the Zakharov Hamiltonian for the gravity water waves and we 
write a truncated approximation of it in action-angle coordinates. As mentionned above, 
we exemplify the technics on the Zufiria 3-modes interaction model. 



2 THE ZAKHAROV k ZUFIRIA HAMILTONIANS 



3 



Section 3 is devoted to the description of the KAM strategy in the general case. We 
identify a part of the Hamiltonian which is small and thus which can be removed by 
a canonical transformation. And the remaining part of the Hamiltonian has a trivial 
trajectory. This method builds the unknown transformation by an implicit equation. This 
change of variables is applied to the set of non-resonant Diophantine initial conditions in 
the neighborhood of a point in the phase space. Then we obtain a trajectory of the original 
perturbed Hamiltonian, in the initial physical coordinates, by computing the inverse of 
this change of variables. 

In section 4 we apply the KAM method to the Zufiria model: we exhibit a part of the 
Hamiltonian which is small for a non-empty set of initial conditions. We also verify that 
the obtained trajectories are not travelling waves (or permanent form) solutions. 
Section 5 is a brief conclusion along with some open questions and work in progress. 
Let us mention that an alternative starting hamiltonian would be any finite-order 
truncation of the ones given in (Dyachenko. 1995), (Craig, 1996) for deep water waves 
or (Craig, 1993) for the finite depth case. We leave this task to further work since it 
would require much more algebra, even though the final result would be physically more 
accurate. 



2 The Zakharov 8^ Zufiria Hamiltonians 

In this section we shall remind some well-known facts of the Hamiltonian formalism of 

surface water waves (Zakharov, 1968). We will study a truncated version of it. 

Let ri{x) be the elevation of the surface wave, ^{x,y) the velocity potential and ■ip{x) 

be the velocity potential evaluated at the surface. The dynamical variables are {r),ip), 

canonically conjugated. We only consider the case where the bottom is at an infinite 

depth. 

Then the Hamiltonian H{ri,ilj) is the total energy of the system given as the sum of the 
kinetic energy and the potential energy Hp, where: 



^6 = 2/0 /-oo y)) + (^2/^(^' y)) 



and 



(2.1) 



1 

Let us first introduce the Fourier representation of r]{x) and ■0(x): 

^(a;) = E%e''" with % = Cfe (2.3) 



i^ix) = ^ke'^'^ with 4 = (2-4) 
feez 

along with the new variables: 

ak^{fik\k\--* +4k\k^)/^/2 (2.5) 
ie 

?7fe = (afe + 61fe)|A;|Vv^ (2.6) 

i^k^{ak-a*_k)\k\~-^/{iV2) (2.7) 
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The dynamical variables are now (17, ip), still canonically conjugated, or also (dfc, d^). The 
Hamiltonian is then reduced to a system of simple equations: 

dH 

dtak = -i^ keZ (2.8) 

By definition, T]{x,t) is a travelling wave if and only if there exists a function g and a 
constant c such that i]{x,t) = g{x — c.t). This is equivalent to: 

dtT] = -c.d^r] or dtfjk = -i.k.c.f]k (2.9) 

or to the requirement that the function 'f)k{t) verifies: 

\/k,t ln{rikit)) = Ck - i.k.c.t, (2.10) 

for some constants e^. 

The expansion of the Hamiltonian in power series of the variables a is formally given by 
(Krasitskii, 1994). 

Many models which are derived from the Zakharov Hamiltonian system have been used 
to study the time evolution of spatially periodic waves, see (Stiasnie & Shemer, 1984), 
(Zufiria, 1987) and (Badulin, Shrira, loualalen & Kharif, 1995). In this paper, we choose 
the 3 waves interaction model introduced by Zufiria, where he assumes that the only 
non-zero variables are 01,02,03 the others being frozen to 0. Due to the existence of an 
extra integral of motion (the total momentum (j2.14j) . besides the Hamiltonian itself), this 
is known to be the minimal model that may exhibit a chaotic behavior. Zufiria gives an 
expansion of the Hamiltonian up to the 4th order of waves amplitudes 0. We can rewrite 
it as follows: 

H = did* + 22d2d2 + 32d3d3 + 2276^did2d3 + d^d2d3^ + 2277|^d^d2 + di^d2) + 

a.a; (ia.a; + 2.,a,a,i + u.a^a^) + a.a^(a,a,i + 2.,,,i^ + ^(a.a^^ + 

278(^did2^d3 + d*d2d3j — 79 (^afat, + d*'^d3^ (2-11) 

/3xi „_9 1 V2 3 + V3 

where 76 = (g j , 7? = 2 * , ^1 = 2 + ' " — 2 — ^ ' 

W3 = 3 + i/-, 78 = 2'i.3-3.(2 + v^+y8) , 79 = ^ (2.13) 
V o o 

Zufiria also notices that this Hamiltonian retains all the resonant terms up to order 4, 
and therefore it can be viewed as the normal form in 3 degrees of freedom. 
For the water waves equation, beside the Hamiltonian, there is an other conserved quantity 
which is the total horizontal momentum given by: 



I = ^ A; dfcd* (2.14) 



fc=i 



(obviously the case 1 = corresponds to a flat surface, for this Zufiria model, where the 
sum is only on positives k). Indeed ()2.14|) is a conserved quantity as it can be checked 
by derivation and the help of the hamiltonian equation ()2.8|) . We can also rederive this 
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fact as follows. Let us introduce the action-angle variables with the following canonical 
transformation from ( I) to {A, e, I, ^, B, defined by: 



0.2 = ^Vl- a-3b- A-3B e^'^ 
v2 

with a > 0, 6 > 0, a + 36 < I 

and A > -a, B > -b, A + 3B<I-a-3b 



(2.15) 
(2.16) 

(2.17) 
(2.18) 
(2.19) 



The cases (a = I, 6 = 0) or (a = 6 = 0) or (a = 0, 6 = 1/3) correspond to the travelling 
waves of class 1, 2 or 3 (respectively) of the Zufiria classification. By studying the stability 
of these solutions, he finds that the most important solution, from a dynamical point of 
view, is the one of class 2. Indeed, for a certain values of the parameter I, the stability of 
this solution changes. 

The change of variables ()2.15ll2.l7jl leads to a new Hamiltonian, still denoted by H: 

H = c + LUiA + LUsB - -fsAB - 74^^ - 75^^ + 



[a + A){I 



3b- A- 3B) 



1 r 

2 



{b + B)2-fQcos{e - (f) + (a + v4) 2 77 cos (26*) 



[a + A){b + B) 



(J-a-3b-A- 35)78 cos(^ + ^) - {a + A)-f9 cos(3^ - ^) 



where: 



c = 

71 

72 

73 

74 

75 
UJi 

UJ3 



7i + I72 + J 

= aui - buj3 + akjs + 0^74 + 6^75 



+ auJi + iu^b 

v2 

V2 + I 



/9 

3(5 + 



2a74 + 673 + ujil + 
a'js + 2675 + UJ3I - 



4V2 



ids 



3 

2' 



1 



1 

V2 



UJ3 



V3 + ^ 



3 

V2 



+ 

(2.20) 

2.21) 
2.22) 
2.23) 

2.24) 

2.25) 

2.26) 

2.27) 
2.28) 

2.29) 
2.30) 



From ()2.2()j) the Hessian of the hamiltonian is non-degenerate in action variables. This 
will be used in section 4 in order to apply KAM method. 

Let us notice that H is independent of ^ and hence we see again that I is a constant of 
motion. 
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3 The KAM method 

We want a canonical transformation which turns an Hamiltonian written in action-angle 
variables, into another one for which we can compute an exact trajectory. And so by 
inversing this transformation we will get an exact trajectory of the initial Hamiltonian 
system. 

So let us consider an Hamiltonian H with L degrees of freedom in action- angle variables: 

(A^)eAxT^ (3.1) 

where T is the torus ]R/27rZ and the domain A is a neighbourhood of a point Aq G M^. Aq 
can be taken to be by translation in the action variable. 

Assume that H is of class with respect to A, and let us write H{A){6) instead of 
H{A,9): this means that H{A) is a function of 9. The Taylor expansion of H around 
A = is then: 

H{A){e) ^c + u.A- f{e) - v{e).A - q{A){e) (3.2) 

where c is some scalar constant, iv is some vectorial constant e R^. Here /, V and q{A) 
are functions of 6 with values respectively in R, R-^, R. 

We denote by ' the derivative with respect to A, and hj d = dg. So that we have 
q{0) = q'{0) = (q is of order 2 in A). We assume that !D°[g"(0)] is an invertible matrix 
from R^ into itself, where D^{f) is the average of a function / over the angles, and: 

D* = 1 - D° (3.3) 

with 1 being the identity matrix. 

If / = \^ = 0, the Hamiltonian equations of motions are: 

dtA = -dH = -dq{A) (3.4) 
dtO^H' ^oj + q'(A) (3.5) 

When A is taken to be zero, we get a trivial equation. And the flow is: 
for any initial condition 9. 

This was the "integrable case", where f — V — 0. Now let us consider the "quasi- 

integrable case", i.e. we assume that / and V are small quantities and we define: 
e = max(|/|, \ V\) for some norm. So the hamiltonian we investigate is quasi-integrable 
near A = 0. To avoid ambiguities we make the hypothesis: 

D0(/) = d0(\/) = (3.7) 

Indeed the quantities D'^{f) and D°(V) can be added to c and u; respectively. 
We want a diffeomorphism T from A x to some domain in R^ x which transforms 
the Hamiltonian H in another one for which we can exhibit some explicit trajectories. 
Let us define the translation operator 7 acting on the action variable by: 

yA e R^ VZ e R^ [7{Z)H]{A) = H{A + Z) (3.8) 
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still valid when the members of this equality are functions of 6. Let us also define the 
dilation operator D (acting on the action variable) by: 



VA G wee Li) 



K^)H]{A) = H{QA) 



(3.9) 



so that ro(e)-i = ©(e-1). 

Let us choose the following ansatz for T: 



T = T2.T1 with 

Ti = T(d - Z) 

T2 = [1 + Gd)-^]D){1 + G) 



and where: 



a : 



a 



da, G = dG 



(3.10) 
(3.11) 
(3.12) 



(3.13) 



so that Z is a constant vector. Let us emphasize that we denote by a dot the derivative 
with respect to 9. Of course we assume that G is small enough to have 1 + G invertible. 
Notice that Ti,T2 (and so T ) are canonical transformations. Indeed: 

TiA = A + a- Z ^ d{TiA) = dA + ad9 (3.14) 
Tie = e =^ d{Ti9) A d{TiA) = dO A {dA + adO) = dO A dA (3.15) 

since d9 A d9 = 0. And instead to prove that T2 is canonical, let us rather prove this for 
Ts = T^^ = D(l + G)-\l + Gd): 

T3A = (1 + G)-^A =^ diT^A) = (1 + G)-^dA + (. . .)de (3.16) 

T^e = (1 + Gd)e = e + G ^ dir^e) = de{i + g) (3.17) 
diT^e) A d{T^A) = de{i + G) A (1 + G)-^dA + de{i + g)a{.. .)de = deAdA 

Q.E.D. 

Notice also that: 



T 



-1 



A 



[1 + G)-\A + Z 

e + G 



a] 



The transformed Hamiltonian H = TH: 

H{A) = (1 + Gd)-^H[{1 + G)A + a - Z] 
can be written similarly as ()3.2|) : 

H(A) = c + iu.A- f - V.A - q{A) 



where c = c — u.Z — T) 



V.a + q{a — Z) 



V.G + (^{a- Z).{l + G) 



u 

/ = (1 + Gd)-^V* [f - uj.a + V.a - V.Z + q{a - Z) 
V = (l + Gd)-^V*\v.{l + G)-uj.G + (^{a- Z).{l + G) 



(3.18) 

(3.19) 

(3.20) 
(3.21) 

(3.22) 

(3.23) 

(3.24) 
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and 

q{A) = {1 + Gd)-^'D* 
Obviously: 



q(a-Z + (l + G).A) - q{a - Z) - q'{a -Z).{1 + G).A (3.25) 



m = g'(0) = (3.26) 



Hence if we manage to find a, G, Z such that / = = then we will get a trajectory 
of the transformed hamiltonian H since the initial data A — remains forever frozen to 
under the evolution of the special hamiltonian: 



H{A){9) = c + u;.A-q{A){9) 



for which the flow is: 




9 





9 + cut 



(3.27) 
(3.28) 



for any initial condition 9. 

Using the inverse transformation T~^, we get a trajectory of the initial Hamiltonian H: 







T 



-1 





9 + ut 



(3.29) 



where: 



T- 



[l + G)-\Z-a) 
9 + G 



(3.30) 



In order to have f — V — we require a, G k. Z io verify: 



uj.G = 'D' 



f + V.a-V.Z + q{a- Z) 
V.{l + G) + q'{a- Z).{l + G) 



(3.31) 
(3.32) 



Let us define the formal inverse of the operator cu.d by: 

r = {uo.dyW* (3.33) 
which is defined on any function of the angles, by its action on the exponential function: 

Vi/ e \ {0} r.e^ = (i.u;.u)-\e^ where e^(9) = e'''-^ (3.34) 

and r.l = 0. 

Of course we need to assume an arithmetic condition on a; namely that it satisfies a 
Diophantine condition: 

Vi/eZ^\{0} \u;.u\>-f\u\-^ (3.35) 

for some numbers 7>0, r > L — 1. The set of such frequencies has a positive Lebesgue 
measure (if 7 is small enough), as is well known. 
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Then ()3.3H) . ()3.32|) can be rewritten formally as: 



a = r 
G = r 

G 



f + V.a- V.Z + q{a - Z) 

V.{l + G)+(^{a- Z).{l + G) 
1 



le: 



l-V[V + q'{a- Z)]d 



T[V + q'{a- Z)] 



(3.36) 
(3.37) 
(3.38) 



Notice that a is defined by a fixed-point equation ()3.36p a = eY{a). We will adjust the 
(vector) parameter Z by requiring that u = u, i.e.: 



V.G + q'{a- Z).{l + G) 







At order e, the equation for Z is: 



rj)0 



q"ma-Z) 



and so: 



D»[g"(0)] D"[g"(0)r/]+O(e 



(3.39) 



(3.40) 



(3.41) 



To solve ()3.36|) & ()3.39|) . we can use a particular version of KAM method (Vittot, 1994). 
See also (Gallavotti, 1994). 



4 The Zufiria model as a quasi- integrable 
Hamiltonian 

Here we apply the previous KAM method to the Zufiria Hamiltonian (j2.2Up . We consider 
I as a parameter (since I is constant on any trajectory). So H becomes an Hamiltonian 
in 2 degrees of freedoms (A, 6*, i?, v?), ie L = 2. We exclude the cases I = 0, or a = 0, or 
6 = 0, or a + 36 = I. So we can write H under the form ()3.2p with: 



f{9,^) = 2 a{I-a-3b) 



6276 cos(^ — ip) + 20277 cos(26') 



+ 



(a6): 



(I — a — 36)78 cos(6' + ip) — 079 cos(36' — ip) 



(4.1) 



Likewise V = V = (Vi, V3) with: 

Vi{e,if) = [b/a] ^ (^-2Q-36) ^^ ^^^^^ -if) + [b/a] ^ (I - 3a - 36)73 cos(^ + <f) 
[I-a-36]5 

. ,si . . N (21 -3a -66) , 
3(a6)579cos(3e - ^) + ^ r^77Cos(2^) 



[I-a-36]: 



(4.2) 



and: 



1^3 (^, = [a/b] 5 cos(^ -f) + [a/b] ^ (I - a - 96)73 cos(^ + if) 



[I-a-36]2 



a[a/6] 279 cos(36' — ip) 



3a 



[I- a -36] 2 



-77 cos(26') 



(4.3) 
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We choose a and b in the domain B defined by: 



{(a, 6) G 



if ; 



a + 36<I; \f\k\V\<eo} 



for the norm: 



I/I 



snp\fii^)\.F^iRM) 
with F<,(z/) = E ^ 



(4.4) 

(4.5) 
(4.6) 



using the Fourier coefficients / of /. Likewise \V\ = \Vi \ + iVsj. 

The parameters a, R, Eq can be taken to be a = 3 (since we have 2 degrees of freedom), 
R = 1/3 and eo = 1/10. B is non-empty when I < 1/40. The shape of this domain is 
shown in the figure below. 

The generators a, G = (Gi.Gs), Z = {Zi, Z^) are given by dHSHI), (ESU), (ESHl)- And 
the solution (|3.29p . (j3.30|) is written as: 





9 

V V J 



I (1 + G) 



V 



-1 
6- 



^3 
6*3 



"3 



(4.7) 



with di = dea , = d^a. Let us recall that the arguments of the functions 

^1} Gri, G3 

are 6, (f. Then the fiow acts on the angle variables by (j3.29p : 

e{t) = 6(0) + ioit (4.8) 
ip{t) = ip{0) + U3t (4.9) 

since a; = by the choice of Z. To the first order in e, we have (cf ()3.36|) . ()3.38|) ): 

a = Tf (4.10) 
Gi = TiVi + Tdef) (4.11) 
Gs = nV, + TdJ) (4.12) 

Now we define the Cantor set K of Diophantine points (a, b) by the conditions (j3.35|) 
where u is actually an affine function of a and b (I being fixed): cf (j2.27|) . (|2.28p . It is 
non-empty and has a positive Lebesgue measure (if 7 is small enough), as is well known. 
Let us also define B = B fl K. We can also prove that it is non-empty and has a positive 
Lebesgue measure (since B is an open set). 

For any value of (a, 6) G B we can write a trajectory of the Hamiltonian (j2.20|) under the 
form P^ - ipn^ & (jTiT|) . And so we get a trajectory of the Hamihonian (jTTT|) : 



ai(t) = 

^3it) = \lb 
1 



Zi - Vdef + 0(£2) exp?(0 + F^i + V^def + i + 0{e^ 



Z3 - Td^f + 0(£2) expz + F^s + V^d^f + 3e + 0(£^ 



02 (t) 



72 



I — aia\ — 35,303 exp(22^) 



(4.13) 
(4.14) 

(4.15) 
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with 



(4.16) 



We have seen in ()2.10j) that ri{x, t) is a travelhng wave if and only if: 

VA;,t hiifikit)) = Ck - i.k.c.t (4.17) 

for some constants 6^, which is equivalent, due to (|2.5|) . to: 

VA;, t ln(afc(t)) = 4 - ^-^-ct (4.18) 

for some other constants Cfc. 

In particular taking k = l^'m order to get a travelling wave, the following should be true: 

^ ln(a + Zi - rSe/) + i{e + T^i + T^def + = ei - ^ct (4.19) 

But this cannot be true since the l.h.s formula is not affine in t because / and V are 
quasi-periodic (cf (j4.ip - (j4.3j) ). Hence our solutions are periodic in x, quasi-periodic in t, 
and not travelling waves (since the Diophantine vector uj is not rational). 



5 Conclusion 

In this paper we considered the problem of existence of solutions of a free surface in gravity 
waves, quasi-periodic in time and periodic in space. We applied the KAM method to the 
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particular Zufiria model of 3 waves interaction which is a truncation of the Zakharov 
Hamiltonian. From the expressions ()4.13|) - ()4.15|) we remark that these solutions are not 
the ones obtained by the Hopf bifurcation, as in (Van der Meer, 1985). For the Hopf 
bifurcation, the solution is obtained by a continuous variation of the parameter, which is 
the celerity, of a permanent form solution. This is not the case for the solutions given 
by the KAM method, they exist in a domain of the phase space which contains the non- 
empty Cantor set B. Moreover, these solutions are not null solutions since a & G are not 
identically zero. And they are not permanent form solutions neither because a and G are 
not constant functions. 

This method may be extended to a finite trucation of the Zakharov Hamiltonian with 
more modes involved. 

In extension to this work, we would like to give a better threshold for the domain of the 
existence of the solution. Our next aim, unsolved for the moment, would be the case of 
a system of infinitely many degrees of freedom. It is known to be a difficult problem in 
KAM theory. 
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